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1 ,_! ' Abstract. - In pump-probe experiments an intensive laser pulse creates non- equilibrium excited 

O ■ electron distributions in the first few hundred femtoseconds after the pulse. The influence of 

non-equilibrium electron distributions caused by a pump laser on the apparent X-ray magnetic 
circular dichroism (XMCD) signal of Ni is investigated theoretically here for the first time, con- 
. sidering electron distributions immediately after the pulse as well as thermalized ones, that are 

not in equilibrium with the lattice or spin systems. The XMCD signal is shown not to be simply 
| ; , proportional to the spin momentum in these situations. The computed spectra are compared 

to recent pump-probe XMCD experiments on Ni. We find that the majority of experimentally 
observed features considered to be a proof of ultrafast spin momentum transfer to the lattice can 
I . alternatively be attributed to non-equilibrium electron distributions. Furthermore, we find the 

XMCD sum rules for the atomic spin and orbital magnetic moment to remain valid, even for the 
Ch ' laser induced non-equilibrium electron distributions. 

o : 
o . 




■ £^ ' Introduction. — Ultrafast control of the magnetiza- 
, ^^tion of ferromagnetic materials is essential for modern 

■ f Technologies such as high-speed magnetic recording and 

. spin electronics [1-3] . Magneto-optical pump-probe tech- 
niqucs have become exploited to investigate the relaxation 
, process and dynamics of the magnetization in the sub- 

■ picosecond regime [4-6] . In such experiments an intensive, 
, • • short laser pulse is applied to optically excite a material 
' . £h subsequently the magnetic response of the system is 

■ probed by a magneto-optical technique, for example, the 

■ Kerr effect. Earlier pump-probe experiments on Ni in the 
'. .^visible range showed a strong reduction of the Kerr sig- 
nal [4], yet, the mechanisms of the magnetization decay 
were not well understood. The possible dissipation chan- 
nels, the role of the spin-orbit interaction (SOI), and the 
influence of the non-equilibrium electron population have 
been discussed [4,7-10]. 

Stamm et al. [11] have recently extended the pump- 
probe technique by using the x-ray magnetic circular 
dichroism (XMCD) to probe the magnetic response of a 
Ni film that was excited through an intensive pump laser 
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pulse. The observed XMCD value at the Ni L3 edge was 
found to be strongly decreased, by 65% of its original 
value, which was interpreted as a quenching of the Ni mag- 
netization within 120 fs. As the XMCD technique gives 
access — through the XMCD sum rules — to separate values 
for the spin and orbital magnetism, a new window on their 
individual response is offered. Stamm et al. [11] concluded 
from a sum rule application that it is the spin compo- 
nent that becomes strongly reduced, whereas the orbital 
moment is practically unchanged, and deduce accordingly 
that the spin becomes reduced within 120 fs through spin 
angular momentum transfer to the lattice. 

However, there are several open questions. The sharp 
drop of the Kerr signal observed in pump-probe experi- 
ments on Ni has been interpreted as evidence for demag- 
netization occurring through SOI within 20 fs [7], but it 
has also been attributed to state-blocking effects [8] , or to 
a combination of these [12]. The state-blocking effect oc- 
curs when certain unoccupied states that are accessible to 
the probe laser under equilibrium conditions are no longer 
accessible, as they have been populated already through 
the pump laser, and also, some of the originally forbid- 
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den transitions become allowed. Magneto-optical pump- 
probe experiments on Ni films [8] indicated that state- 
blocking effects had disappeared from the transient Kerr 
response at about 250 fs after the beginning of the pump 
pulse. More advanced magneto-optical pump-probe po- 
larimctry experiments performed on CoPt3 films demon- 
strated that state-blocking effects had vanished from the 
magneto-optical response 150 fs after the onset of the 
pump [13,14]. The somewhat different relaxation times 
for CoPt 3 and Ni could be possible, because CoPt 3 ex- 
hibits a much stronger spin-orbit interaction than Ni and 
a very different band structure. Therefore, the relaxations 
in these materials have to be considered separately. 

Here we study computationally the influence of non- 
equilibrium effects on x-ray absorption and XMCD spec- 
tra. We show that such effects do considerably modify 
both the x-ray absorption spectrum (XAS) and XMCD 
spectrum, each in a specific way. The predicted modifica- 
tion of these spectra for the case of Ni is found to be in 
good agreement with observed experimental features [11]. 
We consider two forms of non-equilibrium electron distri- 
butions, first, specific modifications of the population of 
states due to the pump laser, which are expected to exist 
to about 100 fs after the pulse, and second, thermalized 
electron distributions, which are expected at about 300- 
400 fs after the pump pulse [15,16]. In the former case, the 
electron distribution is not a Fermi-Dirac (FD) function, 
while in the latter case, the thermalized electrons can be 
described by a FD distribution, but the electrons are not 
yet in equilibrium with the lattice. 

Computational methodology. In our computa- 
tional approach we use, first, experimental information 
as well as ab initio calculated transition probabilities to 
determine the amount of population or depopulation of 
Ni's spin-polarized band states by the pump laser. These 
modified occupation numbers are subsequently used in rel- 
ativistic calculations of the XAS and XMCD spectrum, 
in which the exchange-splitting of the core levels is in- 
cluded. The XAS and XMCD asymmetry are determined 
from the dielectric tensor whose elements we compute 
from the Kubo linear-response expression. This compu- 
tational scheme, combined with accurate relativistic band 
structure calculations, in which the exchange-splitting of 
the core states is taken into account, provides theoretical 
XAS and XMCD spectra that agree very well with exper- 
iments [17, 18]. Our evaluations of the XAS and XMCD 
have been implemented on the full-potential, relativistic 
WIEN2k method [19], employing the local spin-density ap- 
proximation (LSDA). The full Brillouin zone in the cal- 
culation comprised 10 000 k points, the fee Ni lattice pa- 
rameter was 3.42 A. Lorentzian broadening of 1.2 cV for 
the L 2 edge and 0.7 eV for the L 3 edge was applied to the 
calculated spectra to represent the core-hole lifetime and 
smoothen small features of spectra below the experimen- 
tal resolution. We note that certain features of the elec- 
tronic structure of Ni are not correctly described within 



the LSDA, in particular the width of the occupied d band 
and the appearance of a photoemission satellite peak at 
6 eV binding energy [20,21]. These discrepancies can be 
corrected when on-site correlation effects are taken into 
account [22]. 

The number of pump-laser excited electrons per atom 
of the studied sample follows from the effective excitation 
density f exc . Its maximum value fl x ^ x reached due to 
the pulse can be estimated from the information about 
the laser and the absorption in the material, using the 
formula [8]: f™ x = f abs (dm) Hid^ > wnere d m is the 
thickness of the thin Ni film, f a b s is the fraction of incident 
light absorbed in the film, Vm the unit cell volume, E p the 
pulse energy, u> the frequency of the pulse laser, and S is 
the area of the spot irradiated by the pump laser. For the 
experimental values given by Stamm et al. we obtain the 
maximum excitation density f^ x = 0.76, where the fact 
that 15% of each laser pulse excites the sample is taken 
into account, dm=15 nm, E p =2 mJ, / a 6s=0.1, huj=l.59 
eV, and 5'=0.75 mm 2 . The fine accuracy of f^. x esti- 
mate might be limited, but what is important here is that 
it determines roughly the range of values for f exc to be 
expected in the system. Our results are provided as a 
function of f exc for a wide range of f exc starting from 0. 
We note that the present value of f^ x is rather high; the 
energy of the employed pump laser E p (2 mJ) is consid- 
erably larger than, for example, that of E p =1.6 nJ which 
was used in an earlier pump-probe study of Ni [8]. To 
avoid very high occupations of bands above Ep we have 
restricted our calculations to f exc values up to 0.6. 

Next, we compute the occupation numbers of relevant 
bands in Ni, following the approach outlined earlier [23]. 
The simplest picture of a non-FD distribution due to a 
pulse with fixed energy Huj is the following: occupied bands 
down to fku below Ep are partially depleted and previ- 
ously unoccupied bands up to fuo above E-p are partially 
populated. The occupation numbers of those bands that 
are accessible to the pump laser are computed from ab 
initio calculated optical intensities. Note that the op- 
tical excitation conserves spin, thus we examine major- 
ity and minority spin transitions separately. The spin- 
resolved occupation of bands f^ a above the Fermi level 
is computed from f^ g = "^"^ fexc, where u>b i(J is the 
spin-dependent weight of transitions to band b, w& the 
probability of transition with spin a, and n(, )CT is the to- 
tal number of spin electrons that fits in the band in the 
considered energy window^. Below Ep there are several 
flat bands with similar transition strengths and a simpler 
approach is adopted here, neglecting the small differences 
between band-resolved transition strengths, but ensuring 

Our ab initio calculations give that close to 100% of all majority- 
spin transitions go into valence band no. 11, while the minority spin 
transitions are divided between bands no. 10, 11, and 12, with re- 
spective weights «Jfc i 52%, 3%, and 44%. For ri(, j(T , the number of 
electrons that fit in the relevant bands, we obtain 0.39, 0.022, and 
0.28 for spin-down electrons in bands no. 10, 11, and 12, respectively, 
and 0.22 spin-up electrons in band 11. 
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that spin-polarization of electrons removed from occupied 
bands is correctly accounted. 

After the initial laser excitation, fast electron-electron 
processes cause electron equilibration on a time scale of 
about 350 fs [16], leading to a FD distribution with a well 
defined electron temperature T e . The electron tempera- 
ture can be computed from the energy deposited by the 
laser. The average energy per excited Ni atom follows 
from (£)= J Ef FD (E,T)n(E)dE, where n(E) is the Ni 
density of states (DOS). Using the computed Ni DOS, we 
obtain a maximal electron temperature of 11 300 K, which 
corresponds to the laser excitation with f^ x . A tem- 
poral evolution of the electron population is expected to 
reduce T e already directly after the laser pulse. In our 
calculations we therefore consider electron temperatures 
ranging from 1000 to 11 000 K. These temperatures are 
higher than the Curie temperature of Ni (Tc = 630 K), 
what is possible since the thermalized electrons are not 
yet in equilibrium with the spins nor with the lattice. The 
magnetic moment on Ni, however, computed for such a 
high-temperature FD function, is reduced from the equi- 
librium moment, because there exists a high spin-majority 
3c?-DOS peak just below Ep whose contribution becomes 
diminished through the FD function. 

Calculated results. — Fig. QJtop) shows the com- 
puted equilibrium XAS spectrum at the Ni L 3 -edge as 
well as the computed non-equilibrium XAS spectra due 
to state-blocking and thermalization effects. Both state- 
blocking and thermalization cause a shift of the calculated 
Z/3 XAS peak to lower energies and also small broaden- 
ing of the peak width. This is the very same modifica- 
tion of the XAS as the one that has been observed ex- 
perimentally [11]. The XAS difference signal of the non- 
equilibrium and equilibrium spectra exhibits a 5 eV wide, 
oscillator-type structure centered around 852 cV which 
is very similar to the one observed experimentally [11], 
Fig. njbottom) . The modification of the XAS spectrum 
can be understood from the depletion and repopulation: 
first, down to 1.59 eV below Ep the band states are de- 
pleted by the pump laser pulse. These states can be ac- 
cessed in the subsequent x-ray transition, which gives rise 
to extra XAS intensity at lower energies. Repopulation 
removes furthermore accessible states above Ep, leading 
to a reduced intensity at higher energies. Due to lifetime 
broadening the effects are distributed over a several eV 
wide spectral range. 

Fig. [5] shows the full XMCD spectrum computed with 
state-blocking or thermalization effects included. State- 
blocking for f exc = 0.6 reduces strongly the height of the 
L3 and L2 peaks and increases the XMCD signal at the 
low-energy sides of both L3 and L2 peaks. Thermalization 
also reduces the XMCD peak at both L edges, but does 
not create extra signal at the low-energy sides. Note that 
in both cases the change of the XMCD peak position is 
almost negligible when compared to the change in the XAS 
peak position (Fig.[T]). 
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Figure 1: (Colour on-line) Top: The XAS I/3-edge spectra com- 
puted for Ni in equilibrium (black line) and non-equilibrium 
XAS spectra based on the theoretical model with excitation 
density f exc — 0.6 (dashed red curve) and assuming electron 
thermalization with T e — 7000 K (dashed green curve). The 
corresponding difference spectra between the two latter and 
the equilibrium spectrum are shown in the middle panel. Bot- 
tom: Difference between experimental L3 edge XAS spectrum 
of Ni 200 fs after the onset of the pump laser and the XAS for 
equilibrium [11]. 



To provide a detailed analysis, we study the response 
at particular energies, denoted A and B in Fig. [TJ which 
correspond to the XMCD maximum (B, at 852.15 eV) 
and to the XAS value equal to one half of the XAS max- 
imum (A, at 851.52 cV in our model) in accordance with 
the experiment [11]. For later discussion we introduce a 
third energy C given by the maximum of the XAS differ- 
ence curve (at 851.09 eV). The experiment [11] revealed 
a strong reduction of the XMCD signal at energy B. The 
computed changes of XAS and XMCD at energy positions 
A, B, and C depending on the excitation density f exc and 
temperature T e are shown in Fig. [31 The XAS at A is grow- 
ing with the excitation density f exc , while the XMCD at 
B is decreasing, in agreement with the experiment. The 
change of these properties is roughly linear in the exci- 
tation density (examined up to the value f exc = 0.6), as 
was suggested earlier [8] . In contrast to the state-blocking, 
electron thermalization leads to an XMCD reduction over 
the whole edge. This happens because electron equilibra- 
tion involves both a magnetization reduction and some 
blocking of optical transitions due to the broad FD distri- 
bution. Noteworthy, it is possible to find such f exc and T e 
that the calculated relative XMCD reduction at B matches 
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Figure 2: (Colour on-line) Calculated XMCD spectra of Ni. 
Black line: equilibrium XMCD, red dashed line, with state- 
blocking for an excitation density f exc = 0.6, and green dashed 
line, with thermalized electrons of T e = 7000 K. 

the experimental value. 

A complete XMCD spectrum allows to obtain the total 
spin and orbital momentum contributions to the magne- 
tization via the XMCD sum rules [24,25]. The LSDA 
spin and orbital moments calculated for Ni in the ground 
state are tos=0.62/ib and m£=0.043 /ib- Applying the 
sum rules (with the usual neglect of the magnetic dipolc 
term), we obtain tos=0.46^b and m£=0.042/iB. Thus, 
the orbital moment is properly given, but the accuracy 
for the Ni 3d spin moment is weaker; an error of up to 
36% with respect to the exact value was reported previ- 
ously [26] . Applying the sum rules to the electron distribu- 
tion for f exc = 0.6, we obtain that ms changes from 0.46 
to 0.44 [i B and m L from 0.042 to 0.037 fj, B - This finding 
agrees (within the sum rules' accuracy) with the fact that 
the electron's spin is preserved in the pump laser excita- 
tion, but the orbital moment shrinks because transitions 
from d to p states are prevalent. The XMCD sum rules 
are therefore valid for the studied excited system with the 
same accuracy as for the ground state, but one must take 
care when the perturbation of the XMCD spectrum is not 
homogeneously distributed over the edges. Although the 
XMCD signal at the L3 peak is reduced 1.6 times in our 
simulation, the area of the whole L3 edge is reduced only 
5%, because of the compensating XMCD signal rise near 
energy C. Consequently, taking the XMCD peak at B as a 
measure for the reduction of the spin moment would un- 
justly suggest a significant demagnetization, which would 
not be observed when the whole spectral area is consid- 
ered. 

Also in the case of electron thcrmalization, we find-from 
an analogous comparison-that the sum rules remain valid. 
Replacing the spectral integral of the sum rules with the 
XMCD value at B is in this case a better approximation, 
as the XMCD signal over the whole edge is reduced. How- 
ever the reduction is not homogeneous. If we consider, for 
example, a T e of 5000 K, a reduction of the XMCD at B 
of 42% is computed, but the reduction at C is only 12% 
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Figure 3: Change of the XAS and XMCD signals (normalized 
with respect to its value without laser pulse) at the energies 
A, B (and C for XMCD) as function of the excitation density 
feme (top) and electron temperature T e (bottom). 

(Fig. [2]). This difference is due to a partial blocking of 
states that results from the high-T e FD function. A direct 
calculation of the moment in the thermalized state gives 
ms=0.39 hb (?7j_l=0.024 fig) which implies a reduction 
of the spin moment by 37%. An application of the in- 
tegral sum rules gives a reduction by 33%, in reasonable 
agreement. 

Discussion of the temporal evolution. In the 

case of state-blocking a linear increase of the XAS at A 
with f exc is predicted, which is due to the shift of the L3 
peak to lower energy. A smaller increase of the XAS at A 
follows from electron thcrmalization with a FD function 
having a high T e (see Fig. [TJ . Conversely, the XMCD sig- 
nal at its maximum (point B) is reduced for both cases of 
state-blocking and high electron temperature, but for the 
XMCD peak value the latter effect is slightly stronger, i.e. 
the XMCD peak maximum is reduced less by the state- 
blocking effect than by the corresponding high-T e state 
(Fig. This observation allows us to predict that the 
crossover from the non-equilibrium, state-blocking state 
to the electron-thermalized state should be accompanied 
by a reverse of the initially increasing trend of the XAS (at 
A) and, contrary to it, an almost unchanged or slowly de- 
creasing trend of the XMCD maximum at B. At larger 
times a slow increase of the XMCD at B is, however, 
expected, due to the electron-lattice equilibration, which 
acts as a competing factor. 

The experiment [11] shows a linear rise of the XAS sig- 
nal at A with time up to about 220 fs after the onset of the 
70 fs pump pulse, after which the XAS starts to decrease. 
A different temporal behavior is observed for the XMCD 
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at B: the XMCD signal reduces linearly with time up to 
220 fs after the onset of the pulse, after which it remains 
roughly constant. This indeed confirms the behavior we 
predicted above for the crossover to electron thermaliza- 
tion. The here presented theory asserts therefore that for 
Ni in the studied experiment [11], state-blocking is preva- 
lent up to about 150 fs after the 70 fs pulse, when the 
crossover to electron thcrmalization has occurred, which 
enforces a reduction of the XAS at A. We note that the 
here deduced relaxation time of 220 fs for the crossover 
to electron thermalization in Ni is consistent with earlier 
pump-probe magneto-optical measurements [8], whereas 
a moderately smaller relaxation time of 150 fs after the 
onset of the pulse was measured on CoPt3 films [13]. 

Conclusions. — In Ref. [11] important conclusions re- 
garding the spin and orbital moment evolution were drawn 
through application of the XMCD sum rules. It is essential 
to note, however, that the change of the spectral integral 
over the whole L3 edge was estimated from the XMCD 
value at only one energy (at point B). From the presented 
calculations we conclude that the change of the spectral 
integral over the whole L3 edge cannot be estimated from 
the XMCD value at only one energy (B) if state-blocking 
is present. It is therefore not possible to decide about the 
role of transfer of spin angular momentum to the lattice 
from the present experimental data [11], missing a com- 
plete femtosecond-resolved XMCD spectrum. 

As we have found that the XMCD sum rules (in their 
integral form) are valid even in the laser excited state, 
a complete XMCD spectrum of laser-pumped Ni should 
thus provide reliable information on the evolution of spin 
and orbital moments. We demonstrate that state-blocking 
and high electron temperature effects lead to the follow- 
ing changes in the spectra: the shift of the XAS to lower 
energy, the shape of the XAS difference spectrum, the rise 
of the XAS at A, and the strong reduction of the XMCD 
peak at B in the first 220 fs after the onset of the pulse. 
This implies that the presence of non-equilibrium electron 
distributions can largely explain the experimental obser- 
vations [11] without the suggested transfer of spin angular 
momentum to the lattice. Our calculations also predict 
the approximate time evolution of these changes related 
to the crossover to thermalization, which we estimate to 
occur at about 150 fs after ending of the pump pulse. We 
propose that other experimental data, especially those ob- 
tained for femtosecond-resolved XMCD at energies near 
position C, together with the presented theory, would al- 
low to separate unambiguously the various contributions 
to the apparent XMCD and pave the way for a precise 
understanding of the nature of the femtosecond demag- 
netization as well as reveal the timescales crucial for the 
demagnetization process. 
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